Enhanced spin Hall effect by resonant skew scattering in orbital-selective Kondo effect 
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The enhanced spin Hall effect in An metal due to the resonant skew scattering is studied with 
first-principles band structure calculations. Especially the gigantic spin Hall angle 75 = 0.1 observed 
recently (T.Seki et al., Nature Materials 7, 125 (2008)) is attributed to the orbital-dependent Kondo 
effect of Fe in the Au host metal, where the t2a-orbitals are in the mixed-valence region while Cg- 
orbitals are in the Kondo limit. The enhanced spin-orbit interaction by the electron correlation in 
the t29-orbitals leads to the gigantic spin Hall effect. Impurities with 5d orbitals are also discussed. 

PACS numbers: 71.70.Ej, 75.30.Kz, 75.80.+q, 77.80.-e 



Spin Hall effect (SHE) is an effect where the transverse 
spin current is produced by the electric field or electric 
current 0, i,!, % % 0, % i, 0, [Til, [H, [H, Q [Hi . It 
does not require any magnetic field or magnets, and is an 
interesting and promising phenomenon for application to 
spin injection and manipulation in spintronics. SHE is es- 
pecially robust and large in metals due to the large num- 
ber of carriers and Fermi energy l^, U, 13, 3 14, IBl- 
For example, SHE in Pt metal even at room temperature 
14 1 is more than 2 orders of magnitude larger than 
that of GaAs Even larger SHE has been recently 
reported in Au/FePt System [l3|, where the spin Hall 



conductivity is 
75 is as large as 



10^ VI- 
^ 0.1. 



cm and the spin Hall angle 



Naively speaking, SHE is the two copies of anomalous 
Hall effect (AHE) for up and down spins, respectively, 
and the knowledge on the latter can be directly trans- 
ferred to the former. When the longitudinal conductivity 
axx is larger than ~ 10^ Q,~^cm~'^ , and the Hall conduc- 
tivity (Th is much larger than /{ha) = 10^ f2~^cm~^ 
with h and a being the Planck constant and the lattice 
constant, respectively, the dominant contribution to the 
AHE is the extrinsic skew scattering [lil [itI. [isl. [l9l. [20t . 
In this case, au is proportional to axx, and hence the 
Hall angle 7 ~ oujoxx is the well-defined measure for 
the magnitude of AHE independent of the impurity con- 
centration. In other words, the Hall angle can be deter- 
mined by examining the single scattering event due to the 
impurity. The typical value of 7 is estimated as the ratio 
~ A/si? with A being the spin-orbit interaction (SOI) and 
Ef the Fermi energy. Usually, ei? is at least a few eV, 
the A is of the order of 10-20 meV for the 3(i orbitals, for 
example. Therefore, the Hall angle is roughly estimated 
as 7 ~ 10~^. This Hall angle 7 can be enhanced by 
the resonant skew scattering by the magnetic impurity 
17 , 1^ ■ Using the Anderson Hamiltonian describing the 
virtual bound state causing the resonant scattering, the 
anomalous Hall angle can be expressed in terms of the 
hybridization energy A and SOI A, and the phase shift 8\ 
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FIG. 1: (color) The skew scattering due to the spin-orbit in- 
teraction of the scatterer and the spin unpolarized electron 
beam with wavevector k with the angle 6 with the spin po- 
larization S{0)n with n = (k X k')/\k x k'\. 



due to the p-wave scattering as 7 ~ (A/ A) sin^i, which 
can be of the order of 10~^ [l7|. Therefore, the spin Hall 
angle 75 = 0.1 is a surprisingly large value, which needs 
to be understood for the design of the gigantic SHE. 

In this paper, we propose that the local electron corre- 
lation and spin fluctuation further enhance SHE com- 
pared with AHE by the explicit first-principles band 
structure calculation with the Au metal as the host where 
the Fermi energy is at the 6s bands. Namely, SHE is not 
the simple two copies of AHE, and the comparison be- 
tween these two offers a unique opportunity to study the 
many-body effect. 

Let us consider the scattering of a spin = 1/2 particle 
by a potential with SOI. The amplitude of the scattered 
wave is given by [lit 



/t(^) = hm T> + *e'^/2(0)| i) 



(1) 



for incoming up-spin and down-spin electrons, respec- 
tively, where 9 is the angle between the wavevectors of 
incident (fc) and scattered ik') waves, and fi{d))[f2{d)) 
corresponds to the spin non-flip (spin flip) scattering am- 
plitudes. The skewness of the scattering is represented 
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by the function 



S{0) 



2lm[f*{9)f2m 
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(2) 



while the function 1(6) = \fi{e)\^ + f2[e)? determines 
the scattering intensity. The meaning of the skewness 
function S{9) is shown in Fig.l. When the unpolarized 
electrons are incident, the scattered electrons are spin 
polarized as S{9)n with n = {k x k')/\k x k'\. By solving 
the Boltzmann equation with this skew scattering, the 
spin Hall angle 75 is given by Q 



is 



J dm {6)3 (0) sin. 
J dni{0){l - cosO)' 



(3) 



where / d^l is the integral over the solid angle. The nu- 
merator in eq.(3) represents the transverse spin current 
produced by the scattering, i.e., the velocity perpendic- 
ular both to k and n with the spin polarized along n, 
while the denominator corresponds to the usual trans- 
port scattering rate. Without the resonance effect, the 
typical value of 7 is of the order of 10""^, i.e., much smaller 
than unity. Actually, the obtained value of 7 in rcf.[^ is 
1/500 for n-type GaAs assuming the screened Coulomb 
potential. The partial wave analysis gives the expression 
for M0) and /2(0) 



MO) = E 
hie) = E 



P/(cos( 
2ik 
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dcos 



:P/(COS0). 



Putting cq.(4) into eq.(3), wc obtain 

')( 



Im 



7s = 7: 



-2iSi 



2 9 sin^ 6^+4 sin^ 5^ -h 3 [l - cos 2 {S+ - 6^)] 

(5) 

where we have assumed that the resonant channel is the 
d-wave (/ — 2) which is subject to SOI and the scattering 
is characterized by the two phase shifts ~ <^,/=2±i/25 
while that {61) for the p-wavc scattering is assumed to 
be spin-independent. Assuming that 5i for the non- 
resonant J?- wave scattering is small ( \Si\ = 0.1 [13 )j 



75 = — 3(5i(cos2(5^ - 
3[l-cos2(<52+-<52-)]; 



cos2(5^)/(9sin^ J+ -I- 4:Sm'^ S~ 
Therefore, the most important 
factor is cos 25^ — cos 2S2 , which is related to the dif- 
ference in the occupation numbers of the J = 2 ± 1/2 
imp urity states induced by SOI through the Friedel sum 
rule 
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Therefore, the local density of states (DOS) for 
the d-electrons determines the magnitude of SHE, which 
can be studied by the first-principles band structure cal- 
culation as described below. 

We take Au as the host metal, where the gigantic SHE 
has been observed [l3| . The electron configuration of Au 
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FIG. 2: (color) Total, site-, orbital- and spin-decomposed 
DOSs of (a) bulk fee Au, 3.1 % Au-vaeancies, (b) 3.1 % Pt- 
ipipurities, (c) 3.1 % Fe-impurities in nonmagnetic state and 
(d) in ferromagnetie state. 

(4) 

atom is 5d^'^6s^, and the Fermi energy is at the 6s6p- 
bands with relatively small SOI. Therefore, it is difficult 
to explain the giant SHE in terms of the intrinsic mech- 
anism, in sharp contrast to the case of Pt with 5d^6s^ 
14 1 . Because in [l^] the spin current is supplied from 
FePt to Au, it is natural to consider the three possibil- 
ities of the imperfections, i.e., (i) Au vacancies, (ii) Pt 
impurities, and (iii) Fe impurities. From the above con- 
sideration, there are several requirements to obtain the 
gigantic SHE; (I) there should be the resonance at the 
Fermi energy due to either the Kondo peak or the mixed 
valance, (II) the orbital angular momentum should not 
be quenched, and (III) the peak must be split due to SOI 
by the energy comparable or larger than the width of the 
peak. We have studied the impurity states in Au host 
metal for these three cases in terms of the local density 
approximation (LDAn27[ plus on-site Coulomb interac- 
tion U (LBA+U) = 5, J = 0.9eV). We used 
the highly accurate full-potential augmented plane wave 
method, as implemented in the WIEN2k code [2^. The 
atomic positions are relaxed in the presence of the im- 
purity. The numerical results presented below have been 
tested for convergence with respect to the energy cut-off 
for augmented plane waves, fc-points used for irreducible 
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TABLE I: Down-spin occupation numbers of the 3d-orbitals 
of the Fe impurity in Au from LDA+U calculations without 
SOI and with SOI. The calculated magnetic moments are: 
mf = 3.39 hb and m*°* = 3.32 hb without SOI, as well as 
mf*^ — 3.19 fiB, rrta'^ = 1.54 fiB and m^"* = 3.27 /is with 
SOI. The mufRn-tin sphere radius Rmt ~ 2.65ao (flo is Bohr 
radius) is used for both Fe and Au atoms. 
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Brillouin zone and supercell size. Figure 2 shows the den- 
sity of states (DOS) in the presence of one impurity in 
the 2x2x2 supercell. Basically the structure extend- 
ing from -8 eV to -2 eV is due to the 5d bands of Au, 
while 6s6p bands of Au are extended all through in this 
energy range with smaller DOS. The change in the DOS 
is almost confined in the range of the 5d bands in the 
case of (i) vacancy (Fig. 2(a)), and is at around -1.5 eV 
in the case of (ii) Pt (Fig. 2(b)). We could not obtain a 
magnetic state in both of these cases even in the LDA+U 
calculation. This conclusion is consistent with the ear- 
lier study by photoelectron spectroscopy [2^. In sharp 
contrast, we obtained the magnetic state as the ground 
state with both the LDA and LDA+U calculations in the 
case of Fe. Accordingly, the DOS for Fe (Fig. 2(d)) shows 
the spin splitting, and the down spin DOS has the sharp 
peak close to the Fermi energy. (We put the DOS for 
nonmagnetic state in Fig. 2(c) for reference.) This means 
that there occurs the valence fluctuation of Fe ion be- 
tween and c?^. This is reasonable since the Fe in Au is 
known as a Kondo impurity but the Kondo temperature 
is as low as 0.4 K. From these results, we conclude that 
only (iii) Fe impurities could the origin of the large SHE 
satisfying the criteria (I). (II), (III) discussed above, and 
we now analyze Fe impurity in Au in more detail below. 

The band structure calculation for the state with lo- 
cal magnetic and orbital orderings corresponds to the 
mean field theory of the Anderson Hamiltonian [2l|. In 
the limit of isolated impurity atom, this mean field the- 
ory gives the correct energy positions of the peaks in 
the spectral function of the single-particle Green's func- 
tion. Of course, there is no symmetry breaking by the 
local electron correlation, and the true ground state is 
the quantum mechanical superposition of the degener- 
ate symmetry breaking states, i.e., quantum fluctuation 
of spins and orbitals occur. Early theories of Kondo ef- 
fect of Fe in Au have estimated the relevant quantities 
a.s 1.4eV, U ^ 5.4eV, J ^ 0.9eV [H, where A is 
the energy broadening of the virtual bound state due to 
the hybridization with the conduction bands, U is the 



on-site Coulomb interaction, and J is the Hund's cou- 
pling energy between different orbitals. The crystal field 
splitting is considered to be small (~ O.leV), and SOI 
A = SOmeV is even smaller. The resistivity measurement 
at room temperature shows a systematic change as the 
valence of the impurity changes as Ti, V, Cr, Mn, Fe, 
Co, Ni, and shows the dip at Mn, while the maximum at 
Fe (Fig. 17 of Ref. [il). This strongly suggests the peak 
in the DOS for Fe, i.e., the mixed valence case, although 
the low temperature properties have been analyzed by 
the Kondo model [2J| . The magnetic susceptibility mea- 
surement shows the S = 2 magnetic moment, while the 
Moessbauer experiment concluded rather 5 = 3/2 [ist . 
These somewhat confusing situation can be resolved by 
the LDOS in Fig. 2(c-d) together with Table I. Table 

1 shows the spin/orbital magnetic moments and the oc- 
cupation number of each d-orbital in terms of LDA -|- 
SOI + U calculation. Note that the obtained values de- 
pend slightly on the muffin-tin sphere radius Rmt and 
U value, but the semi-quantitative conclusion does not 
change. First, the nonmagnetic state in the LDA calcula- 
tion shows almost no crystal field splitting (of the order of 
O.lcV), which is consistent with the earlier result. How- 
ever, the inclusion of U changes the situation dramati- 
cally, and the eg — t2g splitting is enhanced to be around 

2 eV as shown in Fig. 2(d). This corresponds to the 
orbital polarization due to the electron correlation, i.e., 
the local version of the orbital ordering not by the crystal 
field but by the electron correlation [26[. Therefore we 
conclude that the orbital-dependent Kondo effect occurs 
for Fe in Au; the eg orbitals are in the Kondo limit, while 
t2g-orbitals are in the mixed valence region. Therefore, at 
temperatures above Tk = 0.4 K, the <2g-orbitals. within 
which the orbital angular momentum is not quenched, 
play the major role in the transport properties, while the 
eg-orbitals determine the low temperature Kondo effect. 

Now we consider the orbital polarization within the t2g 
states due to the SOI. Note that i2g-orbitals behaves like 
Ics. = 1 states with xy,zx + iyz, and zx — iyz orbitals cor- 
responding to m = 0, m = 1, and m = —1 states, respec- 
tively. Therefore, SOI is effective within the t2g states 
leading to the energy splitting between the Jeff. = 3/2 
and Jeff. = 1/2 states. Naively, the orbital polarization 
is determined by the competition between the hybridiza- 
tion energy A and the energy splitting due to SOI. Here, 
one must carefully distinguish between the many-body 
states and the single-particle states. The energy sepa- 
ration between the many-body state with different total 
angular momentum J is typically the order of SOI, and 
much smaller than U or J. However, once the many-body 
ground state of d-electrons is fixed, the separation of the 
single-particle state energy, which is the energy difference 
between the TV-electron and N ± 1-electron state, can be 
as large as U or J, which is important for the conduction 
electrons which comes in or out to the c?-orbitals. There- 
fore, it is possible that the electron correlation U plays an 
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essential role, and even a SOI much smaller than the hy- 
bridization energy can produce the large orbital magnetic 
moment = 1.54^^*^ as shown in Table I. Correspond- 
ingly, the single-particle m = 1 state is almost occupied 
(0.819) while m = — 1 state is almost empty (0.087), as 
is seen also from the inset of Fig. 2(d). Figure 2(d) and 
Table I also explain why the AHE is rather small com- 
pared with SHE. As for AHE, the phase shift for to = 1 
(jTi = —1) are almost tt (0), and both of them does con- 
tribute a little to the scattering, while to = state is 
almost in the unitary limit with 7r/2-phase shift. There- 
fore, we do not expect the enhanced AHE, leading to the 
small anomalous Hall angle compared with that for SHE 
17|. For the SHE, on the other hand, such a cancela- 
tion does not occur. One needs to treat the spin/orbital 
fluctuation in this case, but a rough estimate for the spin 
Hall angle can be obtained as follows. For the conduction 
electrons, the energy difference between Joff = 3/2 and 
JcB = 1/2 is that of TO = 1 and to = — 1 states in the 
mean field theory when one considers the Ising type cou- 
pling IzSz- This energy is larger than the hybridization 
energy and we expect the difference of the phase shifts 
between these two channels of the order of tt, giving a 
large spin Hall angle of the order oi Si =0.1 as observed 
experimentally [l3|. This picture for the Kondo effect 
in Fe impurity in Au is different from the conventional 
one [2^ . Further experimental studies using the spectro- 
scopies such as STS and ARPES are highly desired to 
clarify the nature of this fundamental problem. 

SHE associated with Kondo effect is expected for other 
ions such as rare-earth impurity [sot . From this view- 
point, we have also studied the case of 5d impurities, i.e., 
Ta, W, Re, Os, Ir doped into Au. In these cases, the 
appropriate U value is around 0.5 eV, and there occurs 
no magnetic state. However, SOI is almost as large as U, 
and hence the large SHE is expected. The results for the 
DOS and the associated phase shifts are shown in Fig. 3 
for W (Fig. 3 (a)) and Os (Fig. 3(b)) as the two represen- 
tative examples. Even without the magnetic moment, 
SHE is enhanced due to the resonant state at Fermi en- 
ergy, and the corresponding spin Hall angle is estimated 
as 7 = 0.2 sin 2^1 where 61 is the phase shift for the p- 
wave scattering, and we expect |7s| = 0.04 in these cases 
in ser ni- q uant it at i vely agreement with the early experi- 
ment fl7| . 

In conclusion, we have studied the enhanced SHE due 
to the resonant skew scattering. Compared with AHE, 
SHE, where the time-reversal symmetry is not broken, 
gives a unique opportunity to study the electronic states 
near the impurity without disturbing its magnetic be- 
havior. We have presented the formula for the spin Hall 
angle in terms of the phase shifts (eq.(5)), and also the 
first-principles band structure calculation to analyze the 
mechanism of the gigantic SHE. We have shown that SHE 
shed a new light on the Kondo effect, which plays the key 
role to enhance SHE and a new picture for the Fe impu- 
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FIG. 3: (color) Site-decomposed DOSs and phase shifts (S) 
of (a) 3.1 % W and (b) 3.1 % Os in Au. 



rity state in Au has been proposed. This leads to the 
material design of the large SHE even at room tempera- 
ture with the potential application to the spintronics. 

We thank Profs. A. Fert, K.Takanashi, N. Kawakami 
and Y. Otani for fruitful discussions. 
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